Macroscopic tunneling of a membrane in an optomechanical double-well potential 
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The macroscopic tunneling of an optomechanical membrane is considered. A cavity mode which 
couples quadratically to the membranes position can create highly tunable adiabatic double-well 
potentials, which together with the high Q-factors of such membranes render the observation of 
macroscopic tunneling possible. A suitable, pulsed measurement scheme using a linearly coupled 
mode of the cavity for the verification of the effect is studied. 
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Optomechanical systems have seen a recent surge in 
experimental and theoretical interest, culminating in the 
cooling of micromechanical oscillators to within a frac- 
tion of a phonon of their quantum ground state p~H4] and 
the reaching of the strong coupling between cavity field 
and mechanical element [5HI]. Despite these successes, 
few experimental demonstrations of their non-classical 
behavior have been achieved so far. Notable exceptions 
include Ref. pQ, which coupled the mechanical oscillator 
to a Josephson qubit to detect the presence of a single 
mechanical phonon, and Ref. [4], which demonstrated the 
asymmetry between up-converted and down-converted 
photons of a probe laser field, an unambiguous signature 
of the asymmetry between phonon absorption and emis- 
sion. Additional proposals to generate and exploit non- 
classical effects in cavity optomechanics include schemes 
to squeeze a motional quadrature of the oscillator [8]- 
[13], perform quantum state tomography [TJ [14] [15], or 
offer alternative ways to engineer non-classical mechani- 
cal states p~6j H7] including most intriguingly perhaps the 
realization of Schrodinger cat states in truly macroscopic 
systems [18H2Q] . 

This paper extends these considerations by exploring 
the possibility to realize and monitor the quantum tun- 
neling of an optomechanical system operating deep in 
the quantum regime through a classically forbidden po- 
tential barrier. The observation of the tunneling of such 
a truly macroscopic object has not been achieved yet, al- 
though theoretical proposals have been made [21]. We 
find that this can be achieved in a "membrane-in-the- 
middle" (MiM) configuration [22j [23] under conditions 
that are close to being realizable in current state-of-the- 
art experiments. We propose a detection scheme based 
on pulsed optomechanics ideas [15] that permits to moni- 
tor the tunneling dynamics through a series of weak mea- 
surements of the membrane position. 

Our approach relies on adiabatically raising a potential 
barrier, whose parameters can be widely tuned, at the lo- 
cation of a mechanical harmonic oscillator. We show that 
the ground state of the resulting double-well potential 
can exhibit tunneling rates several orders of magnitude 
larger than the decoherence rate of the mechanical mem- 
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FIG. 1. (Color Online) Membrane-in-the-middle optomechan- 
ical system. The dashed lines illustrate the adiabatic double- 
well potential induced by the field mode. The membrane in 
shown in a superposition of left- well state and right- well state. 
The lower figure sketches a scheme for the detection of mem- 
brane tunneling. 



brane, and that a weak optomechanical position mea- 
surement is enough to monitor the tunneling. Besides 
tunneling, the proposed scheme allows for the study of 
the quantum Zeno effect [24] in a mechanical context and 
provides a comparatively simple scheme for the prepa- 
ration and characterization of non-classical mechanical 
states of interest for quantum metrology and sensing. 

We consider a MiM optomechanical system [22 , 23 } [25 ], 
[26] consisting of two fixed mirrors and a partially reflect- 
ing membrane of center-of-mass oscillation frequency ujm 
and effective mass M between them, see Fig.[l] A unique 
feature of this system is the option to realize either linear 
or quadratic optomechanical couplings, depending on the 
precise membrane's equilibrium position. Our proposed 
scheme involves one cavity mode that couples quadrati- 
cally to the membrane to realize a double well potential, 
and a second mode that couples linearly and is excited by 
short optical pulses to first prepare and then monitor the 
membrane position. The mode frequencies as a function 
of membrane displacement x have been verified exper- 
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imentally [6 and resemble very closely sine and cosine 
curves [22j [23] for large enough mode numbers. 

Consider first the quadratic coupling mode. The cor- 
responding system Hamiltonian is 



H = Hi 
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where we have in harmonic oscillator units (ujm — M 
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Here Hm is the energy of the mechanical resonator and 
H op t is the optomechanical interaction, with uj c the res- 
onance frequency of the cavity and g 2 the quadratic op- 
tomechanical coupling constant. The optical pumping at 
frequency u p with pumping rate 77 is described by ^ pump . 
The dissipation terms H K and H 1 account for cavity and 
mechanical damping with rates n and 7 respectively. 

The Heisenberg equations of motion for the membrane 
position can be cast as the second order equation 
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where £ (t) is the noise operator associated with thermal 
damping of the oscillator, and the evolution of the cavity 
field is governed by 
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Given stable parameters, the cavity field approaches its 
steady state in a timescale ft -1 . For n ^> (cjm? 7) it there- 
fore follows adiabatically the position of the membrane 
(note that this neglects the optomechanical membrane 
cooling), and the membrane dynamics are robust against 
fluctuations of the cavity field. Furthermore, a fast cav- 
ity decay rate destroys the quantum correlations between 
the light field and the membrane, in which case it is suf- 
ficient to treat the optical field classically, which we do 
in the following. The intracavity intensity as a function 
of the membrane position is 



I(x,t) 



(^/2) 2 + A 2 (x,t)' 



(7) 



where A (x,i) = S c + g 2 (x 2 ) (t) , with S c = uo c — uo p the 
detuning of the pump from the cavity resonance. In- 
serting this expression into Eq. ([5| and integrating the 
right hand side gives the effective potential acting on the 
membrane, 
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FIG. 2. (Color Online) Ratio of barrier height to ground state 
energy (solid line, left axis) and half minimum separation x m in 
in units of the zero-point width x zpt (dotted line, right axis) as 
functions of the pumping rate 77. Here the coupling strength 
is = — 2 x 10 -4 and the detuning is zero. 



The light field adds an arctangent function to the har- 
monic potential, which for positive coupling g 2 leads to 
tighter confinement and for negative couplings, which 
is the case we are interested in, to a symmetric bar- 
rier whose parameters depend on 77 and S c . For D = 
— n 2 — I6r] 2 g 2 > 0, the potential becomes a symmetric 
double well with two local minima at 
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separated by a barrier of height [27 



= — \x 2 n{n + arctan (2S C I k) + arctan (Vd/ k) 

(9b) 

Typical parameters allow for a wide range of separa- 
tions between minima and barrier heights, as illustrated 
in Fig. [2j which shows x m [ n and as a function of the 
pumping rate 77. The distances between minima can be 
substantial, and the traps should be accordingly shallow 
to give reasonable tunneling rates. For the numerical cal- 
culations, we assumed the quadratic single photon cou- 
pling g 2 /uJM = — 2 x 10 -4 , which for a typical membrane 
with ujm/2tt w 100 kHz translates into a single photon 
optomechanical coupling of 20 Hz. Such relatively high 
coupling rates can be obtained by using avoided cross- 
ings between higher transverse modes [6 . In addition, 
a noteworthy aspect of this scheme is its versatility, as 
weak couplings and/or high decay rates can be compen- 
sated by adjusting the detuning or increasing the input 
power, see Eqs. ([9|, which for the presented calculations 
is of the order of 1 fiW. Generally speaking, low effective 
membrane mass and frequency are desirable to keep the 
tunneling rates high. 

To estimate the tunneling rates in the double-well po- 
tential, we consider the membrane Hamiltonian T~L = 
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YliLi EiC\Ci, with eigenenergies Ei determined the time 
independent Schrodinger equation [p 2 /2-\-U (x)]ip = Eip. 
If Ei and E2 are smaller than the barrier height E^ , their 
values lie very close together. The corresponding eigen- 
states ipi and ip2 are symmetric and antisymmetric, re- 
spectively, but exhibit similar squared amplitude. Thus 
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FIG. 3. Double well and resulting splitting of energy lev- 
els. Also shown the squared amplitudes of the two local- 
ized wavef unctions. Parameters used: (^/cjm = — 2 x 10 -4 , 
77/cjM 176.785. 

the states 

^rA x ) = V^(*W ± ^ 2 (a0) (10) 

are located predominantly in one of the two potential 
wells. A typical situation is depicted in Fig. [3j which 
shows a double-well potential, its energy levels as well 
as the squared amplitudes of the located wave functions 

1pL/R- 

Reexpressing the membrane Hamiltonian in terms of 
the localized modes 

c RtL (x) = l/y/2(c 1 ±c 2 ) (11) 

we obtain 

n = E ^4^+^(4^+4^)+E^- (12) 

je{L,R} " i>2 

with E L = Er = (Ei + E 2 )/2 and the tunneling rate 
J = E2 — Ei . Figure [4] shows the normalized tunneling 
rate J/ojm as a function of the pumping rate t]/ujm, il- 
lustrating its exponential decrease for increasing 77, i.e., 
increasing well separation, see Fig. [2] If more eigenvalues 
Ei lie below the barrier they can be split up analogously 
but we find that ground state tunneling rates become 
unpractically small in such a situation. 

In order for tunneling not to be destroyed by deco- 
herence effects the value of J must be large compared 
to ojm/Q- Recent experiments with SisN4 membranes 
demonstrated Q-factors in excess of 10 6 [28] at frequen- 
cies of 500kHz, with ojm/Q — 0.5 sec -1 and effective 
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FIG. 4. (Color Online) Tunneling rate J of the two lowest 
lying states, normalized to cjm, as a function of the scaled 
pumping rate t]/ujm for — —2 x 10 -4 , k/oom — 10 and 

S c = 0. Note the logarithmic y-axis. 



mass of ~ 50ng. For such a membrane, we find tunneling 
rates in the range of J ~ 100Hz, indicating that the real- 
ization of macroscopic optomechanical tunneling should 
be within reach. 

We now turn to a discussion of the detection and mon- 
itoring of the tunneling. Continuous and precise moni- 
toring of the membrane's position is possible in principle 
through a homodyne measurement of the output light 
from a linearly coupled mode. The effect of such a con- 
tinuous position measurement has been studied in the 
context of single atom tunneling [29 , and reveals two 
main effects. Firstly, continuous measurements can pre- 
vent the tunneling from happening altogether due to the 
quantum Zeno effect. The information about the mem- 
brane's position continuously projects it towards a posi- 
tion eigenstate and thus suppresses the tunneling. Sec- 
ondly, the back-action from the measurement can heat 
up the membrane and thus eventually supply enough en- 
ergy to cross the barrier thermally, without tunneling. 
Although demonstrating the quantum Zeno effect for a 
macroscopic membrane is certainly noteworthy, we fo- 
cus here on monitoring the dynamics of the tunneling 
process. To this end, instead of a continuous measure- 
ment we propose to perform a discrete sequence of weak 
measurements using the pulsed optomechanics technique 
recently introduced by Vanner et al [T5] . 

The quadratically coupled resonator mode is not an 
adequate probe, both because it does not force the mem- 
brane into either well due to the ambiguity of the sign 
of x. Instead, we use a second resonator mode that is 
coupled linearly to the membrane. We assume in the 
following that the center-of-mass motional mode of the 
membrane is initially prepared in its ground state, a step 
that is of course essential and that could be achieved ei- 
ther cryogenically or via an initial optomechanical cool- 
ing stage, depending on the specifics of the experimental 
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setup. Following that preparation, a symmetric barrier 
is raised adiabatically by turning on the power of the 
quadratically coupled mode, so that the membrane finds 
itself in the symmetric ground state of the double well 
potential with density operator p = A first 

measurement pulse then determines the membrane's po- 
sition through a homodyne measurement, its outcome 
x res projecting the membrane with equal probabilities in 
the left or right well. That measurement should be just 
strong enough to project the state of the membrane un- 
ambiguously into one of the two localized states to allow 
for the monitoring of subsequent tunneling. The post- 
measurement density operator after such a measurement 
is 



with 



p(t = 0) 
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where a 2 is the uncertainty of the measurement, inversely 
proportional to the measurement strength [15 . It is de- 
termined by the coupling strength of the light mode to 
the membrane and the pulse intensity. 

Following that initial preparation, the density operator 
evolves unitarily, with the membrane starting to tunnel 
from the random initial well to the opposite one. One 
way to demonstrate tunneling would be to wait a time 
of the order of J -1 before the next measurement pulse, 
assuring the highest probability to find the membrane 
in the opposite well. To obtain information about the 
tunneling dynamics, though, it is preferable to perform 
a sequence of measurements weak enough to keep the 
membrane's energy below the barrier. The outcome of 
a numerical simulation of this process for the situation 
of Fig. [3] is summarized in Fig. [5] For this simulation, 
a sequence of 20 measurements were performed in the 
time J -1 . The insert shows a histogram of outcomes 
of the weak measurements of just one time, indicating 
the spread in values of the membrane position character- 
istic of weak measurements. We note that throughout 
the simulation, the energy remained below the barrier 
height. One can clearly make out the quantum jumps 
in the position of the membrane. If one were to repeat 
the measurement sequence a large number of times, post- 
selecting only those initiated in the same potential well 
would lead to the emergence of the coherent tunneling 
dynamics. 

In summary, we have investigated the possibility to 
observe the tunneling of a macroscopic optomechanical 
membrane in a cavity. The potential barrier is a conse- 
quence of the quadratic optomechanical coupling of the 
membrane to a lossy cavity mode. This set-up is a ver- 
satile tool to create wide, shallow double- well potentials. 
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FIG. 5. (Color Online) Outcomes of a typical sequence of 
20 position measurements on the membrane initially in its 
ground state. The solid line gives the coherent evolution 
(i.e. without further measurements) of (x) for this partic- 
ular trajectory. The histogram depicts the distribution of 200 
measurement outcomes at ujyit = 20, where only trajectories 
which started in the left well were post-selected. Parameters 
of the potential are the same as in Fig. [3] the measurement 
strength is 1/a = 0.02. 



Using parameters from recent experiments, we find that 
achievable tunneling rates can exceed the decoherence 
rates of the membrane by several orders of magnitude. 
The effect could be verified using weak pulses of lin- 
early coupled light. Besides the possibility of verifying 
quantum tunneling at unprecedented scales, it would also 
be a demonstration of the preparation of a non-classical 
"Schrodinger Cat" state. Further studies will include in- 
vestigations of different ways to characterize such me- 
chanical states nondestructively and the potential ex- 
ploitation of their non-classical nature in applications 
such as high-precision quantum metrology, the study 
of decoherence mechanisms, and the quantum-classical 
transition. 
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